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General features of the formalism describing hydrodynamic evolution of 
transversally thermalized matter possibly produced at the very early stages 
of ultra-relativistic heavy-ion collisions are presented. Thermodynamical 
consistency of the model is emphasized. The covariant formulas for the 
moments of the phase-space distribution function are derived. The simple 
model for the transition from purely transverse to standard perfect-fluid 
hydrodynamics is proposed. 
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1. Introduction 

It is commonly accepted that the evolution of the hot and dense matter 
created in heavy-ion collisions at the RHIC energies can be reasonably well 
described in the framework of the relativistic hydrodynamics of a perfect 
fluid [TlElElllElElElElEKin]- Nevertheless, the standard hydrodynamic 
approach, assuming three dimensional (3D) thermalization of the system, 
encounters severe problems concerning initial conditions. The main issue in 
this respect is to explain the very early thermalization time of the created 
system, which is required by the standard hydrodynamic models to describe 
the elliptic flow coefficient and the particle transverse-momentum spectra. 

Recently a possible solution of this puzzle has been proposed by Bialas 
et al. in [TT]. In this approach one assumes that at the very early stages the 
hydrodynamic description applies only to transverse degrees of freedom, 
where the thermalization is quite easy to achieve, while in the longitudi- 
nal direction the motion is reduced to simple parton free streeming. The 
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thermalized, two-dimensional (2D) objects corresponding to the group of 
particles moving with the same rapidity are called (transverse) clusters. 

Originally, the idea of the purely transverse thermalization of matter 
created in heavy- ion collisions was introduced by Heinz and Wong [12\ [TB] , 
with the conclusion that such a model cannot describe the RHIC data. 
However, the new formulation of this concept in Refs. [T3] showed that 
the idea of the purely transverse thermalization may be consistent with the 
experimental results. 

In this paper we concentrate on the discussion of the formal aspects of 
the transverse-hydrodynamics model, which were first studied in [15]. In 
particular, we extent formal results obtained in [111 [H] by including the 
effects of finite parton masses and quantum statistics. 

2. Hydrodynamical equations 

The relativistic hydrodynamic equations of the perfect fluid valid for 
transversally thermalized matter follow from the energy-momentum conser- 
vation law, d^T^" = 0, with the energy-momentum tensor of the form 

r 

Here r = \Jt'^ — is the longitudinal proper time and no describes the 
density of clusters in rapidity. 

Since clusters are 2D objects, whose thermodynamic properties should 
be described by the proper 2D thermodynamic variables, we introduce the 
quantities £2, P^^ ^2 and S2 which denote the 2D energy density, pressure, 
particle density, and entropy density, respectively. The essential feature of 
our treatment is that those variables satisfy the standard thermodynamic 
identities. In our considerations the number of particles is not conserved, 
hence, throughout the paper we restrict ourselves to the case of vanishing 
chemical potential, /x = 0. In this case we have 

£2+P2= TS2 (2) 

and 

de2 = Tds2, dP2 = S2dT, (3) 

where T stands for the temperature. 

The definition of the energy-momentum tensor ([T]) includes also the two 
mutually orthogonal four- vectors and V^. The timelike four- vector U^, 
defined by the equation 



= {uQCOshr],Ux,Uy,Uosmhr]), 



(4) 
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describes the four- velocity of the fluid element and may be obtained from 
the four- vector = [u^, u±,0) = {l,Vx,Vy,0) by performing the Lorentz 
boost along the z axis with the spacetime rapidity 

The quantity is the normalized to unity four-velocity of the fluid el- 
ement in the rest frame of the cluster and = 7 = (l — W'^) ^ is the 
Lorentz gamma factor, determined by the hydrodynamic transverse flow 

The special role played by the longitudinal direction in the transverse 
hydrodynamics, as the direction of the free streaming of the clusters, entails 
the appearance of the additional spacelike four- vector V^, which is defined 
by the equation 

= (sinhr/, 0,0, cosh??). (6) 

In the rest frame of the cluster, where r] = 0, one finds = (0, 0, 0, 1) . The 
straightforward calculation of the components of the energy-momentum ten- 
sor in the rest frame of the fluid element (where in addition to the condition 
r/ = we have also u± = 0) yields the diagonal form 



/ £2 \ 

P2 

P2 

V / 



(7) 



Vanishing of the T"^^ component describing the longitudinal pressure is a di- 
rect consequence of the appearance of the term V^V^ in ([1]) and it indicates 
that there is no momentum transport across the cluster, which also means 
that there is no interaction between the clusters. Furthermore, as compared 
to the standard hydrodynamics, the lack of the longitudinal pressure effects 
in a more rapid transverse expansion of the matter, since all the energy 
accumulated in the clusters at the very beginning of the collision has to be 
used for transverse acceleration only. Of course, it is expected that such a 
situation cannot last long since the collisions between partons tend to ther- 
malize the system also in the longitudinal direction. Hence, within at most 
few fermis after the collision, the transition from transverse to standard 3D 
hydrodynamics is expected. 

One may check by the direct calculation that the energy-momentum 
conservation law leads to the entropy conservation law 



U,^^T^''' = Td^S'' = 0, 



(8) 
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whose structure is similar to that known from the standard hydrodynam- 
ics. In Eq. ([8]) we introduced the entropy current which is given by the 
expression 

= ^saC/'^. (9) 
r 

Equation ([8]) agrees with our expectations, since the evolution in the cluster 
is entropy conserving and the longitudinal expansion of clusters is just free 
streaming, which cannot produce entropy. In other words, the assumption 
that the considered system is a simple superposition of the 2D expanding 
perfect fluids implies directly the adiabaticity of the flow. 

Similarly to the standard hydrodynamics we can incorporate Eq. dH]) into 
the energy-momentum conservation law to derive the relativistic analog of 
the Euler equation 

Wd^iTW") = d'^T + VV^'d^T. (10) 

In comparison with the standard hydrodynamics, an additional term pro- 
portional to V^V^ is present in Eq. (jlOp . By performing projections of 
Eq. ()10p on the four-vectors and we check that Eq. ()10p contains 
only two independent equations. Altogether with Eq. ([8]) we have a system 
of three independent equations for four unknown functions 

S2{t,X±), T{t,X±_), Ux{t,X±), Uy{T,X±). (11) 

Hence, as in the standard hydrodynamics, we need to add the equation of 
state connecting S2 and T to solve the system of hydrodynamic equations. 
The reduction of the number of independent equations from three to two 
is expected, because the evolution along the z axis is fixed and only the 
transverse evolution needs to be determined. 

Our hydrodynamic equations can be further rewritten in a different form 
which is more convenient for numerical calculations. In this case we use 
standard cylindrical coordinates, i.e., we introduce the distance from the 
beam axis r=^r^ -|- r| and the azimuthal angle (p = tan^^ (ry/rx), where 

x± = {rx^Vy). We also introduce the parameterization of the fluid velocity 
in the form 

Vx = v cos(a -|- (/)), Vy = V sm.{a + (j)) . (12) 

The dynamically changing angle a describes deviations of the flow from the 
radial direction. Thus, it specifies the asymmetry of the flow present in the 
most common non-central collisions. 

With the help of the variables introduced above, and using thermo- 
dynamic identities only, we obtain the explicit form of the hydrodynamic 
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equations 

d d d 

— {rs2Uo) + — {rs2Uovcosa) + — {s2Uovsma) = 0, 

d d d 

— (rTuov) + r cos a— {Tuq) + sin a— {Tuq) = 0, 

^ r, f da V sin a \ dT cos a dT , . 

Tul.[- + — ) = (13) 

The quantity d/dr is the total time derivative 

d ^ _|_ ^ _|_ ■^sina 9 ^^^^ 

dr dr dr r dtp 

It turns out that the parameter no, treated as a constant so far, may be an 
arbitrary function of the spacetime rapidity r/. We can show this feature by 
multiplying the original energy-momentum tensor Tj^^^i by uq = (77) and 
incorporating it into the energy-momentum conservation laws 



no{v)T:=, =0. (15) 



Since the tensor Tj^^^i is conserved, Eq. (fT5l) is reduced to the equation 

T„7^i5^no(r?) = 0, (16) 

which, as the straightforward calculation shows, is always fufilled. This 
feature indicates that our model does not have to be necessarily boost in- 
variant, see Fig. [TJ In fact, a closer analysis shows that different initial 
conditions can be applied to different clusters (placed at different values of 
77). This is so because the clusters do not interact and evolve independently 
of each other. 

Although the full 3D thermalization of the system at the very beginning 
of the evolution is hard to achieve, as the system expands the thermalization 
in the longitudinal direction seems to be inevitable. Such situation requires 
the transition from the transverse hydrodynamics to standard hydrodynam- 
ics. Probably, the simplest description of such a transition is achieved with 
the assumption that it has a sudden character. In this case the transition 
may be realized by assuming the so called Landau matching conditions at 
the transition point 

T^^-U, = T^^U,. (17) 

In Eq. (I17p we used the energy-momentum tensor of the perfect fluid valid 
for standard hydrodynamics 

T^^ = {e + P)U^'U^ -Pg^"". (18) 



hydrodynamics 



Fig. 1. Rapidity dependence of the cluster density uq. 



In Eq. (lisp the energy density e and pressure P are 3D thermodynamic 
quantities. Eq. (jl7p . with the appropriate tensors used, may be reduced 
to the local conservation of the energy and momentum (at the transition 
point) 



It must be supplemented by the assumption of the entropy production as- 
sociated with 2D 3D transition, 



where we introduced the 3D entropy density s. 

We note that the 2D 3D transition has been recently addressed in 
the framework of the dissipative hydrodynamics in ITT] . In this frame- 
work the transition is naturally associated with the entropy production. On 
the other hand, the 2D — > 3D transition was also studied in Refs. [TU] 
in the framework of the entropy-conserving anisotropic hydrodynamics and 
magnetohydrodynamics, respectively. Moreover, the Landau matching con- 
ditions similar to those advocated above were also used recently to model 
a transition from the early 3D parton free streaming to 3D perfect-fluid 
hydrodynamics in Ref. [20] . 




(19) 



r 




(20) 



3. Thermodynamics of two-dimenssional systems 



Thermodynamic variables used in the transverse hydrodynamics are 
obtained directly from the potential 0, defined for 2D systems of non- 
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interacting bosons (upper signs) or fermions (lower signs) 

niT, V2,f^)= ±VgTV2 I ^ In (l T e(^-"^)/^) . (21) 
The particles move in the transverse plane whose area is V2 , and their energy 



is £p = m± = ^m? The factor Vg denotes internal degrees of 

freedom. In the case of gluons = 16. Using Eq. ()2ip we can define particle 
density n2 = N2/V2, pressure P2, entropy density S2 = S2/V2, and energy 
density £2 = E2/V2 by using the standard thermodynamic definitions: 

. (22) 



and the relation 

E2 + P2V2 = TS2+fJ.N2. (25) 

Simple analytic formulas may be obtained for massless fermions and 
bosons 

i^gTrT^ 3^gC(3)r3 . 

n2 = ; — , £2 = — (fermions), (26) 

24 47r 

n2 = \^ , £2 = -^^^ — (bosons). (27) 
12 vr 

Here C is the Riemann zeta function. For the two statistics we also find 
that sound velocity Cs is constant and much higher than in the standard 3D 
hydrodynamics, namely 

In the case of the Boltzmann statistics, the analytic expressions can be 
derived also for finite masses 

n2 = ^(m + T)e-™/^, (29) 

ZTT 

P2 = ^(m + T)e-'"/^, (30) 

ZTT 

S2 = ^[m2 + 3mT + 3T2]e-™/^, (31) 
27r 

£2 = ^[r2 + (m + r)2]e— (32) 

ZTT 
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For the classical statistics and massive particles the sound velocity is given 
by the formula 

2 ^ r(m2 + 3mT + 3T^) 



In the limit m Eqs. (|29|1 - (133]) are reduced to those given first in 
Ref. dU. 

4. Phase-space distribution function and its moments 

Our formalism is based on the covariant form of the phase-space dis- 
tribution function that is factorized into the longitudinal and transverse 
part, 

F = f\\g. (34) 

The structure of the transverse part, describing the hydrodynamic evolution 
of a single cluster, comes from the assumption of the local equilibrium 

9iP'U,) = (35) 

where 

p^Ufj, = mj_uo cosh(?/ — r]) — p± ■ u±. (36) 

The longitudinal part follows from the assumption about the free streaming 
of particles 

5{y - ri) no 
f\\ = no = —6 p^V^ , 37 

and realizes the condition y = r] hj using the Dirac delta function. By 
joining together Eqs. (I35p and ()37p we obtain the covariant formula 



F = !^6{p.V)g{p-U). (38) 
r 



Having the explicit form of the phase-space distribution function (I38p at 
our disposal we can calculate its moments, i.e., the particle current N^, the 
energy-momentum tensor T^'^ , and the entropy current S^, defined by the 
momentum integrals 
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The structure of Eqs. ()39p - (I4ip follows from the Lorentz structure of 
the distribution function. Having in mind its scalar character, we find the 
following decompositions in terms of the four- vectors U^^ and V^, 

N^' = aV + bUf", (42) 

T^-^ = a'U^'U" + ^0^"" + c'V^V'' + ^{U^V'' + U''V^'), (43) 

= a"V^' + b"U'', (44) 

where in addition the metric tensor g^^'^ = diag (+1, — 1, — 1, — 1) has been 
used. Since all the coefficients in Eqs. dH]) - (03]) are scalar quantities, they 
can be calculated in the rest frame of the fluid by performing the proper 
projections of Eqs. ()42p - (I44|) on the four- vectors and V^. One can 
check that after simple algebra we obtain Eqs. and ([9]). 

5. Summary 

In this paper we have discussed the general formalism of transverse 
hydrodynamics. The thermodynamic consistency of the model has been 
demonstrated in the straightforward calculations using only thermodynamic 
identities. The hydrodynamical equations have been derived and written in 
the form useful for numerical calculations. The explicit form of the mo- 
ments of the phase-space distribution function has been obtained with the 
method based on the covariant tensor decomposition. The transition to the 
three-dimensional hydrodynamics has been proposed that uses the Landau 
matching conditions. With the appropriate implementation of the transi- 
tion in the numerical code we hope to successfully describe the physical 
observables measured at RHIC. 
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